FiLE  COPY 


MBSIR  76-841 

AN  ANALYTICAL  AND  EXPERIMENTAL  DETERMINATION 

OF  THE  CUTOFF  FREQUENCIES 

OF  HIGHER-ORDER  TE  MODES  IN  A  TEM  CELL 


and 

Myron  L.  Crawford 

Electromagnetics  Division 
Institute  for  Basic  Standards 
National  Bureau  of  Standards 
Boulder,  Colorado  80302 


June  1976 


Prepared  for: 

Electronics  Systems  Division  (AFSC) 
Hanscom  Air  Force  Base,  Mass.  01731 


John  C.  Tippet 
David  C.  Chang 

University  of  Colorado 


NBSIR  76-841 


AN  ANALYTICAL  AND  EXPERIMENTAL  DETERMINATION 

OF  THE  CUTOFF  FREQUENCIES 

OF  HIGHER-ORDER  TE  MODES  IN  A  TEM  CELL 


John  C.  Tippet 
David  C.  Chang 

University  of  Colorado 

and 

Myron  L.  Crawford 

Electromagnetics  Division 
Institute  for  Basic  Standards 
National  Bureau  of  Standards 
Boulder,  Colorado  80302 

June  1976 


Prepared  for: 

Electronics  Systems  Division  (AFSC) 
Hanscom  Air  Force  Base,  Mass.  01731 


U.S.  DEPARTMENT  OF  COMMERCE,    Elliot  L.  Richardson,  Secretary 
Dr.  Betsy  Ancker-Johnson ,  Assistant  Secretary  for  Science  and  Technology 

NATIONAL  BUREAU  OF  STANDARDS,  Ernest  Ambler,  Acting  Director 


a 
* 


CONTENTS 


Page 


FOREWORD  iv 

1.  INTRODUCTION   1 

2.  FORMULATION  OF  THE  INTEGRAL  EQUATION   2 

3.  SMALL  GAP  APPROXIMATION  FOR  THE  GREEN'S  FUNCTIONS   6 

4.  SOLUTION  OF  THE  APPROXIMATE  INTEGRAL  EQUATION   9 

4.1  Cutoff  Frequencies  of  the  Higher-Order  Modes   9 

4.2  Field  Distribution  of  the  Higher  Order  Modes  12 

5.  EXPERIMENTAL  DETERMINATION  OF  THE  RESONANT  FREQUENCIES  OF  TEM  CELLS   13 

6.  RESULTS  AND  CONCLUSIONS  14 

7.  REFERENCES  17 

APPENDIX.      COMPUTER  PROGRAM  FOR  FINDING  THE  ROOTS  OF  AN  EQUATION 

USING  THE  BISECTION  METHOD  18 

LIST  OF  FIGURES  AND  TABLES 

Figure  1.  Cross-section  of  a  TEM  cell   19 

Figure  2.  Nomograph  for  determining  the  cutoff  wavelength  of  the  first 

higher-order  perturbed  mode  in  a  TEM  cell  given  the  ratios  g/a 

and  b/a  20 

Figure  3.  Functional  diagram  of  an  automatic  network  analyzer   (ANA)   system-  21 

Figure  4.  Typical  data  obtained  from  an  automatic  network  analyzer   22 

Figure  5.  Comparison  of  Mittra  and  Itoh's  exact  solution  and  our  approxi- 
mate solution  for  the  cutoff  wavelength  of  the  first  higher- 
order  perturbed  mode  in  a  TEM  cell  23 

Table  1.     Calculated  and  measured  resonant  frequencies  of  the  first  higher- 
order  perturbed  mode  in  nine  different  NBS  TEM  cells  24 


iii 


FOREWORD 


This  report  describes  the  theoretical  and  experimental  determination  of 
the  cutoff  frequency  of  the  first  higher-order  mode  for  several  TEM  cells  6f 
differing  geometry  developed  at  the  National  Bureau  of  Standards    (NBS) .  In 
addition,  the  field  distribution  of  the  first  higher-order  mode  is  found 
explicitly . 

The  effort  is  part  of  a  program  sponsored  by  the  Electronics  Systems 
Division,   Hanscom  Air  Force  Base,   under  contract  number  Y75-917  with  NBS. 
The  purpose  of  this  effort  is  to  evaluate  the  use  of  TEM  cells  for  determining 
the  radiated  emission  characteristics  Of  small  electronic  devices. 

The  theoretical  analyses  were  developed  by  staff  from  the  University  of 
Colorado  under  contract  with  NBS.     The  experimental  work  was  performed  at 
NBS.     Myron  L.  Crawford  of  the  Electromagnetics  Division  was  the  technical 
monitor  for  NBS  and  Charles  E.  Wright  of  the  Electronics  Systems  Division  was 
the  technical  monitor  for  the  Air  Force.     The  period  of  performance  covered  by 
this  report  extends  from  November  1975  to  April  1976. 

Results  of  this  effort  will  be  used  to  develop  higher^order  mode  suppression 
or  rejection  techniques  to  increase  the  useful  frequency  range  bf  a  TEM  cell. 
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AN  ANALYTICAL  AND  EXPERIMENTAL  DETERMINATION  OF  THE  CUTOFF  FREQUENCIES 
OF  HIGHER-ORDER  TE  MODES   IN  A  TEM  CELL 


John  C.   Tippet,  David  C.   Chang  and  Myron  L.  Crawford 

ABSTRACT 

In  addition  to  the  TEM  mode,   higher  order  TE  and  TM  modes  can 
propagate  in  a  TEM  cell.     Since  only  the  TEM  mode  is  desired,  the 
higher-order  modes  restrict  the  useful  frequency  range  of  the  TEM 
cell.     In  this  paper,   the  cutoff  frequency  of  the  first  higher- 
order  mode  is  obtained  both  analytically  and  experimentally  for 
several  TEM  cells  of  differing  geometry.     The  difference  between 
the  experimental  and  theoretical  results  is  shown  to  be  only  a  few 
percent.     In  addition  the  field  distribution  of  the  first  higher- 
order  mode  is  found  explicitly. 

Key  words:     Shielded  strip  line;   TEM  cell. 

1.  INTRODUCTION 

In  our  previous  report  on  the  TEM  cell,  an  effort  was  made  to  investigate 
the  properties  of  the  dominant  TEM  mode   [1] .     In  this  paper,   the  first  few 
nigher-order  TE  modes  in  the  TEM-cell  will  be  studied.     They  are  of  interest 
since  they  determine  the  upper  frequency  limit  below  which  only  the  TEM  mode 
can  propagate.     Thus,  knowledge  of  these  modes  may  enable  us  to  increase  the 
useful  frequency  range  of  a  TEM-cell,   for  example,  by  placing  selectively 
absorbing  material  in  the  cell  at  locations  where  the  effect  is  least  on  the 
propagation  of  the  TEM  mode. 

As  is  obvious  from  the  sketch  shown  in  figure  1,   the  TEM  cell  can  be 
viewed  as  a  special  case  of  a  shielded  strip  line  for  which  the  strip  is 
embedded  in  a  homogeneous  dielectric.     Solutions  leading  to  the  cutoff  fre- 
quencies of  the  higher-order  modes  in  shielded  strip  line  structures  usually 
involve  a  numerical  search  for  the  roots  of  an  infinite  determinant  obtained 
from  a  mode  matching  at  the  interface  between  the  partial  regions  separated 
by  the  center  septum  in  the  cell.     A  problem  arises,   however,   since  the 
infinite  determinant  must  be  truncated  in  order  to  obtain  numerical  results. 
As  Mittra  has  shown    [2],   the  partitioning  of  the  set  of  equations    (i.e.  the 
number  of  modes  used  in  each  of  the  partial  regions)   must  be  chosen  in  the 
correct  ratio  in  order  that  the  solution  of  the  truncated  set  of  equations 
converges  to  the  correct  value.     In  fact,   an  incorrect  choice  of  the  par- 
titioning ratio  leads  to  a  violation  of  the  edge  condition  [3]. 

In  order  to  bypass  the  difficulties  associated  with  relative  convergence, 
an  alternative  approach  was  sought  in  this  report.     An  integral-equation  formu- 
lation for  the  electric  field  in  the  gap  was  found,   and  an  approximate,  but 
analytical  solution  of  the  integral  equation  was  obtained  using  the  assumption 
of  a  small  gap  between  the  center  septum  and  side  walls.      (All  of  the  TEM  cells 
currently  in  use  at  NBS  indeed  have  small  enough  gaps  so  that  the  analysis 


applied.)     This  allows  us  to  represent  the  cutoff  frequency  of  the  first  higher 
order  mode  as  a  perturbation  to  the  TE^q  mode  in  a  rectangular  waveguide.  The 
method  developed  in  this  paper  is  applicable  to  TEM  cells  for  which  the  center 
septum  is  displaced  vertically  from  its  center  position,  and  may  also  apply  to 
TEM  cells  for  which  the  regions  above  and  below  the  center  septum  have  dif- 
ferent dielectric  constants   (i.e.  a  shielded  strip  line).     In  addition  to  the 
cutoff  frequencies,   the  electric  field  distribution  of  the  higher-order  modes 
can  be  found  explicitly. 

In  order  to  check  our  solution,  we.  compared  our  results  to  those  obtained 
by  other  authors   [4,5,6].     In  particular,  Mittra  and  Itoh   [4]  used  a  singular- 
integral  equation  formulation.     Their  results  reduce  exactly  to  ours  when  a  small 
gap  approximation  is  used  in  their  equations.     However,  the  results  obtained  by 
Brackelmann,  et  al.    [5]   and  Bezlyudova   [6]  were  not  useful  for  comparison  since 
their  result  for  a  zero- thickness  septum  reduces  to  the  cutoff  frequency  of  the 
unperturbed  TE^q  mode  in  a  TEM  cell,  independent  of  the  width  of  the  center 
septum.     If  we  regard  the  upper  and  lower  regions  of  the  TEM  cell  as  two  rec- 
tangular waveguides  coupled  through  a  small  aperture,  then  the  use  of  the  so- 
called  coupled  mode  theory  predicts  that  two  system  modes   (symmetric  and  anti- 
symmetric)  can  exist,  both  of  which  have  approximately  the  same  propagation 
constant.     The  symmetric  mode  is  unaffected  by  the  edge  condition  at  the  ends  of 
the  septum,  while  the  anti-symmetric  mode  is  more  strongly  affected  by  the 
width  of  the  septum   [7].     It  is  clear  that  the  results  of  Brackelmann  et  al .  [5] 
and  Bezlyudova   [6]  were  intended  for  the  symmetric  mode.     However,  because  the 
excitation  of  the  TEM  cell  is  basically  of  a  TEM  nature  which  yields  an 
oppositely-directed  electric  field  in  the  upper  and  lower  regions  of  the  cell, 
it  is  more  relevant  to  obtain  results  corresponding  to  the  anti-symmetric  case. 

The  resonant  frequency  of  the  first  higher-order  perturbed  mode  was  also 
obtained  experimentally  in  this  paper.     Using  the  well-known  relationship  between 
the  resonant  frequency  and  the  cutoff  frequency,  excellent  agreement  was  obtained 
between  the  theoretically  and  experimentally  measured  resonant  frequencies. 

2.      FORMULATION  OF  THE  INTEGRAL  EQUATION 

A  cross-section  of  a  TEM  cell  with  an  off-center  septum  is  shown  in 
figure  1.     The  strip,  of  width  2w,   is  located  a  distance  g  from  each  vertical 
sidewall  and  is  assumed  to  have  negligible  thickness.     In  addition  it  is 
embedded  in  a  homogeneous  dielectric  of  permittivity,   eQ  and  permeability,  yQ. 
The  outer  wall  of  the  cell,  of  width  2a  and  height  2b  =  b1  +  b2 ,  and  the  strip 
are  assumed  to  be  perfect  conductors.     The  region  above  the  strip  is  designated 
region  1  while  that  below  the  strip  is  designated  region  2. 
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We  will  attempt  in  this  paper  to  obtain  only  the  cutoff  frequencies  of  the 
TE-type  modes  since  the  TE^g  mode  is  usually  the  lowest  order  mode  next  to  the 
dominant  TEM-mode  in  the  cell.     For  the  TE  modes  it  is  well  known  that  the 
fields  in  region  1  can  be  found  in  terms  of  an  auxiliary  scalar  function, 
¥ ^  ,  which  satisfies  the  following  relations. 

(V2  +  K2)Y(1)   =0  in  S,,  region  1  (1) 


and 


where 


3  =o  on  the  metal  walls,   i.e.,  (2) 


n 


B-C  and  A-F-E-D 


K2  =  k2   _  y2;  (3) 

K  =  -2>Vo  =  (4) 

and  we  have  assumed  that  the  fields  propagate  as  e  1U)t-1Yz>  ^s  tjie  transverse 

Lapiacian,  which,   for  our  geometry  is  given  by  32/3x2  +  82/9y2,  and  3n  repre- 
sents the  outward  normal  derivative.     cj  is  the  angular  frequency,   A  is  the  wave- 
length, k     is  the  wave  number  and  y  1S  the  propagation  constant  which  is  equal 

°  (1) 
to  zero  at  cutoff.     In  addition  to   (2),  we  will  also  require  that  3  ¥        be  non- 

(1) 

zero  in  the  gap  regions,   that  is,  A-B  and  C-D.     If  9  *        were  zero  there  also, 
then  we  could  write  the  solution  to   (1)   immediately  as  just  an  unperturbed 
rectangular  waveguide  mode  unaffected  by  the  presence  of  the  center  septum. 
In 

follows 


In  terms  of  ¥       ,  the  field  components  in  region  1  can  be  written  as 


H     =  K2Y(1)i"  (5) 
z  z 

and 

Et  =  iwyQVt(T(1)  )x  i"z  (6) 

where  a    is  a  unit  vector  in  the  z-direction  and  V     represents  the  transverse 

z     —     3         —  3 
gradient   (a    -r—  +  a    -r— )  . 
3  x  3x        y  3y 

In  order  to  determine  Y ^ ,  we  will  write  it  as  a  superposition  of  a 

complete  set  of  basis  functions  \l>^'   as  follows 

Tmn 

In  S,  =     I     hllUil?  (Xfy)  •  (7) 

1  L      mn  mn 

m,  n 

Ama    are  the  unknown  amplitudes,  and  the  basis  functions,   <p^'  ,  must  satisfy 
the  following 

+  «di>2»ii'J  =  °  (8> 

t        mn  mn 
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and 


9  *(1>  = 
nrmn 


=  0 


(9) 


on  £, ,   the  boundary  of  S,  ,  i.e.  A-B-C-D-E-F.     It  may  appear  that   (7)   is  logically 

(1)         .  '  (1) 

inconsistent  since  the  basis  functions,  satisfy   (9)  while  their  sum,  4* 

rmn 

satisfies  3  f        /  0  in  the  gap  regions.     As  shown  by  Friedman   [8],  however, 

(1)  n 

f        is  discontinuous  at  y  =  0  and  thus  the  limit  y  =  0  and  the  sum  in  (7) 
cannot  be  interchanged.     In  other  words,  the  requirement  of  9  ¥  ^  0  at  the  gap 
should  be  understood  as  applied  to  the  plane  y  =  0+  where  0+  =  e  and  e  ■*  0. 
The  solution  of   (8)   and   (9)   is  obviously 


\b  (1)  (x,y)  = 
rmn  1 


2  1 

h 

(x+a) 

cos 

L2a 

cos 

J 

(10) 


with  given  by 

mn    ^  J 


mn 


©2  + 


(11) 


and  the  normalization  constant  appearing  in   (10)   is  included  for  convenience. 

Using  the  orthogonality  properties  of  the  basis  functions,   the  amplitudes, 

A*     ,  can  be  written  as  follows 
mn 


A(1)    =  A  A  / 
mn  m  n  ^ 


(x-  #y,)*I^)  (x' ,y«)  ds' 


(12) 


where 


A.  = 
l 


i  =  0 
otherwise 


(13) 


In  order  to  find  an  alternative  expression  for  the  surface  integral  appearing 

Subtracting  the  resulting 


in   (12),  we  multiply   (1)   by  and   (8)  by 


equations  and  integrating  over  S^,  we  obtain  the  following 


(K 


(1)  2 
mn 


-  K2)   /  ds  =  /  N 

fal  bl 


(DV2¥(1)  .  ,(1)^(1)]  ds 
mn    t  trmn 


(14) 


Using  the  two-dimensional  form  of  the  divergence  theorem  and  Green's  second 
identity   [9],  the  right-hand-side  of   (14)   can  be  written  as 


R.H.S.  = 


rmn    n  nrmn 


(15) 


The  second  term  in  (15)  is  zero  by  virtue  of  (9).  Thus  the  unknown  amplitudes 
given  in   (12)   reduce  to 


(1)  _ 


A  A 
m  n 


mn 


7     rmn  n 


(16) 


Because  of  (2),  however,  the  line  integral  in  (16)  only  contributes  in  the 
gap.  Thus, 


4  *(1)3 
*     rmn  n 


dl  =  - 


-w  a 

/  +  / 
-a  w 


\b  (1)  (x,0)  3  *  (1)  (x,0)  dx 
ymn  y 


(17) 


since 


3  yd)  =  _  |_  yd)  at  y  =  0 
n  3y  1 


(18) 


In  order  to  combine  the  two  integrations  appearing  in  (17),  we  will  consider 
two  cases:      (1)   the  odd  modes  for  which  ¥(1)(x,y)   =  -4f (1 )  (-x  ,y)  ;  and   (2)  the 


(1) 


(1) 


even  modes  for  which  4*  (x,y)  =  yVJ"'  (-x,y)  .  The  following  analysis  will  apply 
for  the  odd  modes  for  which   (17)   can  be  written  as 


a  ^d)3  yd)  dl  _  _r  3  ^u/  (  0)  [4/-^  (x  o)  -  rii  (-x,0)  ]  dx. 
*       rmn    n  y  1Tmn       '  rmn  ' 

£.  w  2 


(1) 


(1) 


(1) 


(19) 


The  term  in  the  square  brackets  in   (19)   is  zero  if  m  is  even.     Thus   (16)  gives 

"VrP     /     3Y(1)  (x,0)^}  (x,0)  dx       (m=l,3,5,  .  .  .) 

(20) 


A(1>  = 
mn 


(K2  -  K(1)    )  w 
mn 


0 


(m=0,2,4,...) 


Equation   (20)   can  now  be  inserted  into  (7)  which  then  determines  the  scalar 

function,   ¥        in  region  1  in  terms  of  the  unknown  cutoff  wave  number,  K,  and 

an  integration  of  the  unknown  electric  field  in  the  gap  (3  ¥  (x,0)«  Ex(x,0)) 
Thus 


(1) 


(x,y)   =  2  I 


~2  '(1)     2,   *2pk,n<*'y>    ^  V(1)(X''0)*2pll,n(x,'0)   dx ' 

p,n    (K^   -  „)       *     '  w 


2p+l , n 


where  p  =  0,   1,  2, 


(21) 


A  similar  relation  can  be  written  in  region  2  by  replacing  all  subscripts  and 
superscripts  2  by  1  and  by  replacing  y  by  -y.     In  the  gap,  we  require  con- 
tinuity of  Hz  and  E^  which  is  equivalent  to  requiring  the  following 


4-  (1)  (x,0)   =  *  (2)  (x,0) 


(w  <  x  <  a) 


(22) 


and 


3  V  (1)  (x,0)   =  3  *  (2)  (x,0) 

y  y 


(w  <  x  <  a) 


(23) 


When  (22)  and  (23)  are  imposed  on  (21)  and  its  equivalent  expression  in  region 
2,  we  obtain  the  following  integral  equation. 


/     3  4*  (x  1  ,  0)Gq  (x,x'  )dx'   =0       (w  <  x  <  a) 

5 


(24) 


w 


wnere 


G   (x,x')   =     I  I 


n 


j=l  p,n  (K' 


K(j)     2.  *2P+1(X'0)*2P+1(X,'0) 
2p+l,n; 


(25) 


and  we  have  dropped  the  superscript  on  3  ¥(x*,0)   since  it  is  continuous  across 
tne  gap  as  well  as  the  subscript  "n"  on  the  basis  functions  since  they  are 
evaluated  at  y  =  0 .     The  Green's  function  given  in   (25)   is  singular  at  x  =  x'. 
Thus,  as  shown  in   [10,  pp.   197-198]   the  integral  equation   (24)   is  to  be 
interpreted  in  the  principal  value  sense. 

A  similar  expression  can  be  obtained  for  the  even  modes.     The  resulting 
integral  equation  is  identical  to   (24)   if  we  replace  GQ  by  Gg  which  is  given  as 

2  A  A 

G   (x,x')   =     I       I   E-Jlpj          *9_(xf0)ip7    (x',0).  (26) 

e  j=l  p  n   (K2-  K^)2   )       2p  2p 

J       tfl  2p,n 

If  the  septum  is  located  in  the  center  of  the  TEM  cell   (i.e.  b^  =        -  b) 
then  the  summation  on  "j"  in   (25)   and   (26)   can  be  dropped  and  replaced  by  the 
factor  2. 

3.      SMALL  GAP  APPROXIMATION  FOR  THE  GREEN'S  FUNCTIONS 

Tne  integral  equation  given  in   (24)  could  be  solved  numerically  at  this 
point  using,   for  example,   the  moment  method.     In  this  method,  the  unknown  func- 

Using  the 

orthogonality  properties  of  the  basis  functions,  an  M  x  M  matrix  equation 
ensues,   and  the  requirement  that  the  determinant  vanish  then  determines  the 
cutoff  frequencies  of  the  higher-order  modes.     However,   since  a  strictly 
numerical  technique  does  not  furnish  much  insight  as  to  how  the  cutoff  fre- 
quencies change  as  the  gap  varies  and  may  also  lead  to  relative  convergence 
difficulties   [11],  we  will  solve   (24)  by  approximating  the  Green's  function 
using  a  small  gap  assumption,  i.e.  Kg  <<  1.     The  solution  to  the  approximate 
integral  equation  can  then  be  obtained  analytically. 

In  order  to  find  an  approximation  for  Gq,  we  note  that  the  summation  on 


tion,   3  f(x',0),  can  be  approximated  by  an  M-term  Fourier  series 


in   (25)   can  be  written  as  follows 


n=0    [K2  -  K^*  2] 
2p+l,n 


b  .1 

2 

_1 

T\  j 

2p+l 


-  I 


n=l  (n- 


_a(j>2; 

a2p+r 


(27) 


wnere 


l2p+l 


b  . 

TT 


-  (2p+l)2(^) 


(28: 
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The  summation  on  the  right-hand-side  of   (27)   is  given  as  [12] 


a2p+l;       ^a2p+l       Z  2p+l 


so  that   (27)   reduces  to 


An  b* 

J   fM  t—  =   ^-r-t        COt(TT  a^-3,.] 

n=0    [K2  -  k<3>     2]  2 it  a^)  2P+1 

2p+l,n  2p+l 


(30) 


Thus,  GQ  given  in   (25)   can  be  written  as 

2    b2      cotdrai^ ,) 

Go(x,x')   =     I     ^  I   (jf^-  ^2p+l(x'0)*2p+l(x,'0)'  (31> 

:=1  2 tt  p  a2J+1 

Tne  summation  on  "p"  in  (31)  cannot  be  done  in  closed-form;  however,  we  notice 
tnat  for  large  "p" 

cot  (ira^l-  )  r_  i  , 

-2£±L_  *  -     2a     _i   (32) 


aipll  V  (2P+1) 

so  tnat  if  we  add  and  subtract  the  right-hand-side  of  (32)  to  Gq  we  obtain  the 
following 

Go(x,x')   =G^1)(x,x')   +G^2)(x,x')  (33) 

wnere 

m                 2    b2.  r.x      g»,      (x,0)  ip„  +1  (x'  ,0) 

j  =  l  2tt  b.J   p  (2p+l) 


and 

,0.  2  b2 

G^2) (x,x«)   -     J     -1  J 

j=l    2  TT  p 


cot  (Tra^i,  ) 


2p+l'   +  2a 


a2Pll  (2P+1)J 


^2p+1(x,0)^p+1(x',0).  (35) 


-2  (2) 

The  term  in  the  square  brackets  in  (35)  now  decays  rapidly  as  p  .  Hence,  G^ 
approacnes  a  constant  whenever  |  x-x 1  |  ■*  0.  On  the  other  hand,  the  terms  in- 
volved in  decay  only  as  p  and  thus  will  have  a  logarithmic  singularity 
as   |  x-x '|       0.     Therefore  under  the  assumption  of  a  small  gap  the  dominant 

contribution  to  G    will  come  from  G^. 

o  o 

Substituting  for  the  basis  functions  from   (10)  ,   and  making  the  following 
change  of  variables 

t  =  a  -  x  (36) 

and 

f   =  a  -  x'  (37) 
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can  be  written  as 

o 


(1) 


(t,f)   =  I 


p=0  (2p+l) 
which  can  be  summed  as  [13] 


cos 


(2p  +  l)TT  (t  +  f  ) 


2a 


+  cos 


G(1)(t,t')   =  ^jln  tan  —  It+t'l 

n  l4a  '  J 


+  In  tan 


(2p+l)TT  (t-t ' ) 

2a 


^-  |t-f 


(38) 


4a 


(39) 


The  largest  value  of  | t+t 1 |  is  2g  so  that  if  we  assume  that  g/a  <<  1  we  can 
replace  the  tangents  in   (39)  by  their  arguments  and  obtain 


(t,f  j 


In 


4a 


•2_4-  >  2 


(I    «  1] 


(40) 


(2) 


In  order  to  obtain  an  approximate  form  for  G^      we  notice  that  if  we 

restrict  attention  to  only  the  first  few  higher-order  modes  for  which  K  is  not 

too  large,  then  the  sum  in   (35)  converges  very  rapidly,  so  rapidly  in  fact  that 

the  arguments  of  the  cosines  appearing  in  the  basis  functions  in   (35)  remain 

(2 ) 

close  to  zero.     Thus,  G^      is  essentially  independent  of  x  and  x'   for  gK  ■>  0 
and  we  can  write 


(2) 


2     .  b  .  » 

■  I  ?  -1  I 

j=l      a  p=0 


cot(irot2pll) 


2p+l 


2a 


(2p+l) 


(41) 


(2gK  <<  1) 


(42) 


Since  we  have  assumed  a  small  gap,  it  seems  reasonable  to  represent  the  modes 
in  the  TEM  cell  as  perturbations  of  rectangular  waveguide  modes;  that  is,  we 
let 


K2  =  K(S)      +  6(S) 


qr 


(q=l,3,5.  .  .) 


(r=0,l,2 


) 


(s=l,2) 


with 


5  (s) 
qr 


<<  K 


(s) 


(43) 


(44! 


and  K 


(s) 


qr 


cannot  be  so  large  that   (42)   is  violated.     Then  the  main  contribution 

to  the  summation  in   (41)  occurs  when  ai^l,    =  n  or  when  K  -  Ki^'-      .     There  may 

2p+l  2p+l,n  2 

be  more  than  one  term  for  which  this  occurs  since  it  is  possible  to  have  two  or 

more  rectangular  waveguide  modes  with  the  same  cutoff  frequency.     We  will  take, 

for  example^  q  =  1  and  r  =  0  which  corresponds  to  a  perturbation  of  the  TE^q 

mode  in  a  rectangular  waveguide  which  is  the  first  higher-order  mode  that  we 

(s) 


would  expect  to  appear, 


and  only  one  term  in   (41)   is  dominant. 

,(2)  , 


K_l~'   does  not  depend  on  "s 

For  this  case, 


qr 


in  this  case  since  r  =  0, 
(2) 


(G, 


TE 


10 


2 

I 

j  =  l 


cot (bj61Q) 


+  2- 


a6 


10 


*2-  I 

TT    p=l     (2p  +  l) 


can  be  written  as 

TTb  . 


1  -  coth 


(2P+i; 


2a 


(45] 
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Tne  summation  on  "p"  in   (45)   is  very  small  compared  to  the  dominant  p  =  0  term, 
so  we  will  neglect  it,   subject  to  a  consistency  check  when  is  determined. 

In  any  case  the  summation  is  rapidly  convergent  and  can  be  evaluated  numerically 
if  needed.  Thus, 


i  +  cot  Vio  +  cot  VlO 


(46) 


TE. 


10  ""10 
In  the  general  case  for  non-degenerate  odd  modes,  we  obtain  in  a  similar  way 


n.  2       fcot(b  -C  ) 

qr         J        [ 1  ^qr 


+  -    L    jl  -  coth  — (2p+l) 

11  p=0    (2p+l)  2a 


(q-li 


(47) 


wnere 


(s)  = 


qr 


(48) 


The  preceding  analysis  has  dealt  with  the  approximation  of  G  .     A  similar 
analysis  also  applies  for  G  ,   and  we  therefore  only  quote  the  results. 


g]1)  (t,f)    -  J  ln[(J)  |t2-t'2|] 

6  a. 


(49) 


and  in  the  general  case  for  non-degenerate  even  modes  we  obtain 


(Ge  >TE(s, 

qr 


2  'cot(b.E 

y    J   L'<*r 


r  (S) 

^qr 


2 

irq 


+  1 


-    1  -  coth  — J— 

it  p=l  p   [  a  j 

p^q/2 


qtrb  11 


+  — 

IT 


,    COt  i— =  1  i 

2  ■  2  a  - 


(50) 


It  is  interesting  to  note  that  all  of  the  dependence  on  position   (through  the 

variables  t  and  t1)   is  contained  in  G^^  while  all  of  the  dependence  on  the  as 

(    )  '  f  2 ) 

yet  undetermined  K  (through  5       )   is  contained  in  G  . 
2  ^       qr  o,e 

4.      SOLUTION  OF  THE  APPROXIMATE  INTEGRAL  EQUATION 


4 . 1     Cutoff  Frequencies  of  the  Higher-Order  Modes 


Inserting  the  approximate  form  for  Gq  into   (24)  we  obtain 


/  f(f)  ia[(fe) 
o 


t2-t'2 |]  +  G^2) J  df  =  0 


(51) 


where 


f  (f  )    =  3yY(a-t'  ,0)  .  (52) 
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If  we  make  the  following  substitutions 


TTt 

4a '" 

TTt' 

4a  '" 


(53) 
(54) 


Uo  =  &  (55) 


.9 

AQ  =  /     g(t')   dt;  (56) 
o 


h(u')  =  f'(^^-);  (57) 


then   (51)   reduces  to 

J~°  h(u')   In   |u2-  u,2|   du'  =  -        Go^)ho'  (58) 


TT2  „(2) 


o 

But  the  solution  of   (58)   is  known  to  be  [14] 


h(u')   =   ^   (59) 

[u2  -  u'2P 
o 

where  C  is  an  arbitrary  constant  which  is  related  to  AQ  defined  in   (56) . 
Inserting   (59)    into   (58)   we  obtain 

fU°   ^  r-  In   |u2-u'2|   du'   =  —  G(2)A  (60) 

o       [u2  -  u'2P  4a      °  ° 

which  can  be  written  alternatively  as  follows 

fU°   -  r  In   |u-u'  I   du  '  =  —  G(2)A  •  (61) 

r2  12!^  a  O  O 

-u     \m     -  u 1 2 ]  2  4a 

o  o 


If  we  make  the  following  substitutions 


and 


then   (61)   reduces  to 


u  =  u     cos  8  (62) 
o 


u'   =  u     cos  d>  (63) 
o 


C  /     ln[u     Icos  6  -  cos  d>|]d<f>  =  G(   '  A  .  (64) 

1  o'  Y|T4aoo 

o 


However,   since    [10,   p.  171] 


00 

In  I  cos  0  -  cos  <j>  |   =  -ln2  -  2     \     —  cos  mG  cos  m<J>  (65) 


m=l 
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(64)   reduces  to 


TrCln 


°l   =  -JL-  G(2)A  .  (66) 
4a       o  o 


Aq  in  (66)  can  be  found  in  terms  of  C  by  changing  variables  and  substituting 
(59)    into   (56)   with  the  result  that 


A     =  2aC.  (67) 
o 


Thus,    (66)   and   (67)  give 


G(2 


u 

o  1  , 

~        +  —  In 

2  TT 


5~    =  0.  (68) 


^2 
(2 ) 

Substituting  for  GQ  and  uq  from  (46)  and  (55)  respectively  we  obtain  the  fol- 
lowing characteristic  equation  for  the  perturbed  TE,q  mode. 

|[cot  b161Q  +  cot  b251Q]   =  [ln(^|)   -  2].  (69) 

Provided  that  the  remaining  sum  in   (45)   is  indeed  small  under  the  small  gap 
assumption,   this  is  then  the  equation  for  determining  &-^q-     In  particular,  when 
the  gap  ratio   (g/a)   is  much  smaller  than  one,  the  perturbation  solution  can  be 
obtained  simply  from  the  small  argument  expansion  of  the  cotangent  to  yield 

r2  Tib/a 

610   =  -  ~     \    ,8a,  (70) 
blb2  ln(ig> 


and  from   (43)   the  cutoff  wavenumber  is  then  given  by 


K 


TT  Trb/a 

2a       b,b0  ln(— ) 
1  2  "rrg 


(71) 


It  is  of  interest  to  point  out  that  the  form  of  K  as  given  by    (71)    is  essen- 
tially the  same  as  the  one  obtained  in    [15]  which  treats  a  similar  problem  in 
which  coupling  occurs  through  only  one  aperture. 

As  a  trial  solution  for  any  numerical  or  graphical  search  of  the  root  of 

(69)  ,  we  note  that  since  the  value  of  the  cotangent  can  vary  from  00  to  0  when 
t>1^10  increases  from  0  to  tt/2    (assuming  b^  >  b^)  ,   the  expected  solution  there- 
fore lies  in  the  range  where  6^^  varies  from  0  to  Tr/2b,  which  is  the  midpoint 
between  the  unperturbed  TE^q  mode  and  the  unperturbed  TE2Q  mode.     For  g/a 
ratios  from  .01  to   .3  and  a/b  ratios  from  .5  to  1,   a  nomographic  solution  of 

(69)   was  obtained  which  allows  one  to  easily  determine  the  cut-off  wavelength 
graphically.     This  is  included  in  figure  2. 

The  integral  equation  with  Ge  instead  of  Gq  is  identical  if  we  replace  "a" 
by  a/4.     Thus,   for  example,   for  the  perturbed  TE2q  mode  we  obtain 

i[cot  bl62()  +  b262Q]   =  L|°[ln(§f)   -  2]  (72) 

from  which  the  perturbed  cutoff  wavenumber  can  be  obtained  as  above. 
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4 . 2     Field  Distribution  of  the  Higher  Order  Modes 

The  field  distribution  of  the  higher-order  modes  can  be  obtained  by  sub- 
stituting the  aperture  field  found  in   (59)   into   (21) ,  but  first  we  can  evaluate 
advantageously  the  unknown  constant,  C,  in  terms  of  the  "voltage"  drop  across  the 
gap.     Substituting  for  f(t')   from   (52)   into   (56)   and  realizing  from  (6)  that 


E  (x',0) 
d  *(x'  ,0)  =   

Y  iwu 


(73) 


we  obtain  for  A  . 

o 


o 


=  -J—— —  J3    E   (x',0)dx'  =  — 

ltOU      J  X  cou 


(74) 


o  a-g  o 
where  V  is  the  potential  difference  across  the  gap.     Thus,  C  in   (67)   is  given  as 

iV 


C  = 


2acoy 


(75) 


and  the  aperture  field  from  (59)   can  then  be  written  as 

iV       r..2  ..,21-^ 


h(u')  = 


[u2  -  u'2] 

2atou  o 


(76) 


Transforming  the  integration  in  (21)  from  x'  to  u'  via  (37)  and  (54)  and 
inserting  h(u')   from   (76)  we  obtain 


(1) 


;2p+l 


2 

h 

{-} 

7T 

^ab^ 

given 

as  [16] 

\ 

j  o  cos  [  (2p+l) 2u' ]du'  (?7) 

o  /uz  -  u' 2 

o 


jUo  cQ8[(2p+l)2u']du'  =  !j   [2(2p+l)u  ] 


(78; 


where  JQ  is  a  Bessel  function  of  order  zero.  Thus  by  analogy  with  (21)  the 
fields  in  regions  1  and  2  can  be  written  as 

A   (-l)P  j   r  (2p+l)g]  sin  [  (2p+l)^x-]  costj^) 


(3 


(x,y)  =  \-*-\\**-\  I 
ab  . ;  yix>\i  '  p,n 


(K2   -  K<j>  2) 
2p+l,n 


The  summation  on  "n"  in   (79)   can  be  written  as  follows 

1 


A  cos 
n 


y   

n=0    (K2   -  K^i,  2) 
2p+l,n 


b 

2 

TT 

-  I 


n?ry 
cos  |  g— *- 


2a(j)2       n=l   (n2  -  a(j)2) 
a2p+l       n  x   ^r  a2p+l; 


(j  =  1,2; 

(79) 
(80) 
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wnere  ai         was  defined  in  (21 
2p+l 

is  given  as  [12] 


The  summation  on  the  right-hand-side  of  (80) 


cos 


cos 


2p+l' 


2p+l 


2p+l 


J  2p+l 


sin(7ra(^+1) 


(81) 


so  that   (79)   reduces  to 

'2bjW;„  i     -     (-l)PJ^[(2p+l)J|]sint(2p+l)^]  cos 


Y  (D)  (x,y) 


Tra 


iY_|  i 

mu  '  p=0 


2aJ 


From   (28)   and   (43)   we  have  that  a 


(j) 


2p+l 
VlO 


kb  .     -  2p+l 


sin(7ra2pil) 


(82: 


for  r  =  0  so  that  restricting 


attention  to  the  perturbed  TE^g  mode  we  have 
2 


(x,y) 


,-   |  cos[610(y+b.)] 


a6lo'  ^yo 


sm 


sin  6 ,  „b  . 
10  3 


2  a 


2b  . 

TT 


iV 


I 

9=1 


(-DPJ0[  (2p+l)  JflsinI  (2p+l)g] 

I™ 
2p+l 


sinUa2p!l} 


(83) 


Since  the  potential  difference,  V,   is  still  arbitrary,  we  can  normalize  V  in  such 


a  way  that  the  amplitude  of  sin  (•=— )   is  unity  with  the  result  that 


where 


(j) 


b  .6. 


sin  6 ,  nb  .  Q 
 10  j 


I  TTX        ,        j  10 

(x,y)   =  sin | —    +  — '  

2aj         tt        cos  [&1Q  (y+bj  )  ] 


(84) 


.   (-1)        [  (2P+l)If  ]  sin  (  (2P+1)  g  cos  [  (^l)  .ag+1 
q  =      I   .  L  


p=l 


2p+l 


sin(TTa2P!l) 


(j  =  1,2: 


:85) 


In  this  form  the  first  term  in  (84)  can  be  recognized  as  the  dominant  TE^q  mode 
of  a  rectangular  waveguide  and  the  second  term  as  a  small  perturbation.  As 


g/a  ->  0, 


10 


0,   and  the  perturbation  term  in   (84)  vanishes 


5 .      EXPERIMENTAL  DETERMINATION  OF  THE  RESONANT  FREQUENCIES  OF  TEM  CELLS 

Measurements  were  made  using  an  automatic  network  analyzer  to  determine 
the  network  scattering   (S) ,   parameters  for  a  number  of  TEM  cells  as  a  function 
of  frequency.     These  S-parameters  comprise  a  total  characterization  of  the  cell 
as  a  two  terminal  network  from  which  the  cell's  complex  impedance,  insertion 
loss,   and  resonant  frequencies  can  be  determined. 
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The  operational  details  of  the  system  used  to  make  these  measurements  has 
been  published   [17]   and  are  not  repeated  here.     However,   the  test  system's 
functional  diagram  is  shown  in  figure  3.     The  system  has  three  main  sections, 
namely,  source,  measurement,  and  computer.     The  signal  generator  provides 
the  RF  power  over  the  frequency  range  of  interest   (0.1  GHz  and  above)  required 
to  pre-calibrate  the  system  and  test  the  unknown  device  or  the  TEM  cells.  The 
frequency  and  level  of  the  RF  source  is  stabilized,   level  controlled,  and  pre- 
selected by  command  to  the  computer.     The  RF  is  applied  to  the  device  under 
test  via  the  S-parameter  test  set.     Then  depending  upon  how        and        are  set, 

the  parameters  S-^,   S22'   S12'   or  S21  are  measured* 

The  "complex  ratio  detector"  measures  the  amplitude  ratio  and  phase  dif- 
ference between  the  reference  and  test  channels.     This  information  is  digitized 
and  routed  to  the  computer  via  the  instrument  interface. 

The  computer  then  takes  the  S-parameter  data  and  stores  it  either  as  calibra- 
tion data,  if  measuring  standards,  or  raw,  corrected  data  if  measuring  the  TEM 
cell.     This  data  is  then  processed  for  correction  and  routed  to  the  displays 
and  to  other  computer  peripherals . 

Summarizing,   the  measurement  sequence  essentially  is  broken  into  two 
steps.     The  first  is  the  preparation  of  the  system  through  pre-calibration  to 
make  the  desired  measurement.     The  second  is  to  measure,  digitize,  and  store 
the  S-parameter  data  in  the  computer  and  then  convert,  correct,  and  output  the 
data  in  the  desired  form.     Examples  of  data  obtained  for  a  typical  TEM  cell  are 
given  in  figure  4.     Resonant  frequencies  are  indicated  by  the  high  VSWR. 
Note  how  narrow  the  bandwidth  is  for  each  resonance.     This  is  indicative  of  the 
small  RF  losses  within  the  cell's  conductive  walls  and  center  plate  and  is 
an  indication  of  the  cell  Q  factor. 


In  this  section  we  will  discuss  the  result  obtained  for  the  resonant  fre- 
quency of  the  first  higher-order,   perturbed  mode  in  a  TEM  cell  with  a  sym- 
metrically placed  center  septum  via  three  different  techniques;   1)   the  experi- 
mental measurement  as  described  in  the  last  section,  2)   theoretical  calculations 
based  upon   (69),  and  3)   a  numerical  result  based  upon  an  exact  integral  equation 
formulation  as  obtained  from  equation   (58)   in  the  work  of  Mittra  and  Itoh  [4]. 
In  both  of  the  last  two  methods,   the  cutoff  frequency,   f  ,  is  actually  cal- 
culated,  and  this  must  then  be  converted  to  the  resonant  frequency,   f res via 
the  following  well  known  result 


6  . 


RESULTS  AND  CONCLUSIONS 


(86) 


where 


c  =  l//y  e 
Ko  o 


(87) 


and  d  is  the  resonant  length. 
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As  shown  in  the  figure  included  with  table  1,   the  two  ends  of  the  TEM  cell 
are  tapered,   so  that  the  resonant  length  is  not  very  well  defined.     As  a  first 
approximation,  however,  we  will  take  the  average  overall  length  as  the  resonant 
length. 

In  order  to  find  the  zero  of   (69)  ,  a  root-finder  based  upon  the  bisection 
method  was  used.     This  is  contained  in  the  Appendix.     This  same  root-finder 
was  also  used  to  find  the  zeroes  of  Mittra  and  I ton's  equation   (58),  which  for 
a  symmetrically-placed  center  septum  and  a  homogeneous  dielectric  filling 
reduces  to  the  following  for  the  TE  modes 

Y1  (0) 


where 


1  — £         K.J  =0 

kl 


K1  =  a2  =  |(1  +  cos  ™)  ;  (89) 


and 

Y,(0)        S'(0)    +  2d,(0)E1 

— *        =  —x   .  (90) 

k1  k1[S(0)   -  21  ] 

The  symbols  in  the  above  equations  are  defined  in  Mittra  and  Itoh's  paper  and 
in  general  must  be  evaluated  numerically.     However,   for  a  small  gap,   it  can  be 
shown  that  the  solution  of   (88)   is  given  approximately  by 

S (0)   -  21    =  0  (91) 

or  just  the  denominator  of  (90)  equal  to  zero.  Note  that  ■*  0  in  (88)  as 
w  ->  a.  Furthermore  (91)  can  be  shown  to  reduce  to  our  result  given  in  (69) 
for  g/a  ■+  0 . 

Shown  in  figure  5  are  the  calculated  cutoff  wavelengths  computed  via 
Mittra  and  Itoh's  eq.    (58)   and  our  result,   eq .    (69)   for  two  different  TEM  cell 
geometries.     As  expected,   the  two  curves  become  indistinguishable  for  very 
small  gaps.     The  w/a  ratios  of  the  TEM  cells  currently  being  used  at  NBS  are 
all  greater  than  0.72.     Thus  figure  5  shows  that  only  a  few  percent  error 
occurs  even  for  those  cells  with  the  largest  gaps.     It  is  also  interesting  to 
note  that  Mittra  and  Itoh's  result  predicts  a  cutoff  wavelength  given  by  that 
of  the  TE^^  mode  in  the  full  cell  for  septum  widths  up  to  nearly  half  of  the 
width  of  the  cell. 
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Shown  in  table  1  are  the  calculated  and  experimentally  measured  resonant 
frequencies  of  the  first  higher-order,  perturbed  mode  in  nine  different  TEM 
cells  used  at  NBS.     The  column  labeled  "Form  Factor  Group"     refers  to  those 
cells  which  are  scaled  versions  of  one  another.     As  can  be  seen,  the  largest 
percentage  error  between  our  approximate  solution  and  the  measured  value  is 
6.4%  for  w/a  =  .72  or  g/a  =  0.28.     For  smaller  gaps   (w/a  =  .83)   the  percentage 
difference  is  typically  1%. 

Thus  we  are  able  to  predict  the  useful  frequency  range  of  a  TEM  cell  within 
a  few  percent.     Using  the  same  technique,  additional  work  needs  to  be  done  to 
extend  our  solution  to  TEM  cells  which  have  different  dielectrics  above  and 
below  the  center  septum,  i.e.,  a  shielded  strip  line.     In  addition,  if  we  re- 
formulate the  problem  in  terms  of  an  integral  equation  for  the  current  on  the 
septum,   then  we  should  be  able  to  obtain  an  equivalent  form  for  small  center 
septums,   i.e.,  w/a  <<  1.     The  formulation  of  the  problem  for  TM  modes  has  yet 
to  be  accomplished.     When  this  is  completed,  we  should  be  able  to  look  at  the 
TEM  mode  solution   (a  degenerate  TM  mode)   and  obtain  a  formula  for  the  charac- 
teristic impedance  of  a  TEM  cell  with  an  off-center  septum.     Thus,  the  solution 
of  the  conventional  strip-line  problem  can  be  obtained  by  the  same  technique 
described  in  this  report. 
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